Scattering properties of weakly bound dimers of fermionic atoms 
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We consider weakly bound diatomic molecules (dimers) formed in a two-component atomic Fermi 
gas with a large positive scattering length for the interspecies interaction. We develop a theoretical 
approach for calculating atom-dimer and dimer-dimer elastic scattering and for analyzing the in- 
elastic collisional relaxation of the molecules into deep bound states. This approach is based on the 
single-channel zero range approximation, and we find that it is applicable in the vicinity of a wide 
two-body Feshbach resonance. Our results draw prospects for various interesting manipulations of 
weakly bound dimers of fermionic atoms. 
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I. INTRODUCTION 

The studies of degenerate atomic Fermi gases attract a 
great deal of interest as they are directed to reveal novel 
macroscopic quantum states and provide links between 
quantum gases and condensed matter systems. Experi- 
ments with ultracold two-component clouds of fermionic 
atoms widely use Feshbach resonances for the intercom- 
ponent interaction (scattering length a) . This allows one 
to switch the sign and tune the absolute value of a which 
at resonance changes from +oo to — oo. On the positive 
side of the resonance (a > 0), one expects the forma- 
tion of weakly bound diatomic molecules, which represent 
composite bosons and can form a Bose-Einstein conden- 
sate (BEC). On the negative side (a < 0), the inter- 
action between atoms of different components is attrac- 
tive and they should undergo the well-known Bardeen- 
Cooper-Schrieffer (BCS) superfluid pairing at sufficiently 
low temperatures. The BCS transition temperature is 
much lower than the Fermi energy, which makes it diffi- 
cult to achieve this transition. 

Molecular BEC for a > and atomic BCS pairing for 
a < 0, describe the system sufficiently far from the reso- 
nance. In the vicinity of the resonance, where the density 
and the scattering length satisfy the inequality n\a\ 3 > 1, 
the gas enters a strongly interacting regime. This BCS- 
BEC cross-over regime has been discussed in the litera- 
ture in the context of superconductivity [3,I2JMB and for 
superfluidity in two-dimensional 3 He films tjJ.|6J. Atomic 
Fermi gases in the strongly interacting regime are ex- 
pected to have a comparatively high superfluid transition 
temperature 0, IU El EI and are characterized by a 
universal behavior of interactions [l2L fl3L IT3 . Ilj| and a 
universal thermodynamics |16| . They are now actively 
bei ng s tudied in relation to the nature of superfluid pair- 
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Recent investigations of two-component 40 K and 6 Li 



atomic Fermi gases are marked by remarkable achieve- 
ments. Those include the formation 124 . 125. 26, and 
Bose-Einstein condensation [23, HiJ Hfl |U |32| of long- 
lived weakly bound diatomic molecules (dimers) on the 
positive side of a Feshbach resonance (a > 0) , and a BEC- 
type behavior of fermionic atom pairs in the strongly 
interacting regime |33l l34j . An anomalous dependence 
of frequencies and damping rates of quadrupole exci- 
tations on the interaction strength, observed in exper- 
iments [35l |3|| , may be a signature for a transition from 
a superfluid to collisionless regime. Strong evidence for 
the superfluid regime was obtained in the Innsbruck ex- 
periment through the radio-frequency measurement 
of the pairing gap for a strongly interacting Fermi gas of 
6 Li. 

The studies of the strongly interacting regime for the 
BCS-BEC cross-over require the knowledge of many- 
body correlations. In particular, one should reproduce a 
correct equation of state in the limit of BEC of a weakly 
interacting gas of dimers for a > (see 22]) . Hence, 
one should know the interaction between these dimers. 
For a large a they are weakly bound and have a large 
size (~ a) which greatly exceeds the characteristic radius 
of interaction between the atoms. In our previous work 
[8Sl | we have outlined a method for studying the elastic 
interaction between such molecules and their collisional 
relaxation to deep bound states. This method is based on 
the zero-range approximation, and the dimer-dimer scat- 
tering length is found to be 0.6a. The imaginary part 
of the scattering amplitude, originating from the colli- 
sional relaxation, is extremely small. Being in the highest 
ro- vibrational state these diatomic molecules are charac- 
terized by a remarkable collisional stability. The phys- 
ical reason is the Pauli principle in combination with a 
large size of the molecular state (small momenta of bound 
fermionic atoms): collisional relaxation is suppressed as 
it requires at least two identical fermions with small mo- 
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menta to approach each other to a short distance J38| . 

In this paper we present a detailed analysis of elastic 
and inelastic atom-dimer and dimer-dimer interactions 
in the zero-range approximation. The paper is orga- 
nized as follows. In Section II we introduce the zero- 
range approximation in a way it has been done for the 
two-body problem (see, for example, [39l HJj). Section 

III contains a general description of the zero-range ap- 
proximation for the case of three particles. In Section 

IV we review the problem of atom-dimer elastic scatter- 
ing in the zero-range approximation, and in Section V 
present a derivation for the relaxation of the dimers to 
deep bound states in atom-dimer collisions. Sections VI 
and VII contain a generalization of the zero-range ap- 
proximation to the case of four particles. In these Sec- 
tions we present a detailed derivation of the results for 
the elastic dimer-dimer interaction and for the relaxation 
of the dimers to deep bound states in dimer-dimer colli- 
sions. In Section VIII we show that our results can be 
used for weakly bound diatomic molecules obtained in 
two-component atomic Fermi gases by using wide Fesh- 
bach resonances, and in Section IX we conclude. 



II. ZERO-RANGE APPROXIMATION FOR THE 
TWO-BODY PROBLEM 

We first follow well known results and introduce the 
zero-range approximation as this has been done for the 
two-body problem (see We consider elastic pair 

collisions between cold distinguishable atoms interacting 
with each other via a spherically symmetric potential and 
assume that their de Broglie wavelength is much larger 
than the characteristic radius of this potential, R e . In 
other words, we have the condition kR e <C 1, where k 
is the relative wavevector of the atoms. In this case the 
scattering is dominated by the s-wave contribution. The 
behavior of the wavefunction at interatomic distances 
r 3> R e is governed by the scattering length a which is 
related to the scattering phase shift as 6 — — arctan(fca). 
Given the scattering length, the details of the interatomic 
potential at distances r < R e are practically irrelevant for 
scattering parameters and give rise to corrections of the 
order of kR e or smaller for the scattering amplitude |4l| • 

The key idea of the zero-range approximation is to 
solve the equation for the free relative motion of two 
atoms placing the Bethe-Peierls boundary condition on 
the wavefunction at a vanishing r: 
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which can also be rewritten as 



if) oc (1/r — 1/a), r — > 0. 



(1) 



(2) 



One then gets a correct expression for the wavefunction 
at distances r R e . 



For the free relative motion of two colliding atoms, the 
Schrodinger equation takes the form 



(V;! + fc 2 W = 0. 



(3) 



A general solution of this equation for our scattering 
problem, which is valid at any finite r, is given by 



tp = exp(ikr) + W?(r,0), 



(4) 



where G(r, r') = (l/47r|r — r|') exp(ik\r — r'|) is the Green 
function representing the solution of Eq. yj with the 
right hand side equal to <5(r — r'). The s-wave part of the 
incident wave given by the first term on the right hand 
side of Eq. is equal to sin fcr/fcr. As the wavefunction 
if) should satisfy the boundary condition (JTJ at r — > 0, for 
the coefficient h we immediately obtain 



h = - 



1 + ika 



exp(2iS) — 1 
2ik ' 



(5) 



Then Eq. (J3J reproduces the well-known result for the 
s-wave part of the wavefunction at r 3> R e - 
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(6) 



The use of the zero-range approximation is especially im- 
portant for the case of resonance scattering characterized 
by the scattering length \a\ 3> R e . Then, for interparticle 
distances in the interval JJ e « r « 1/fc, from Eq. © one 
finds that Eq. (J2J gives a correct result for if) at distances 
of the order of or even much smaller than \a\. 

For a large positive scattering length a» R e , there is 
a weakly bound state of two atoms. The binding energy 
£o and wavefunction <fio of this state at distances r R e 
can also be found in the zero-range approximation. 

The free relative motion of atoms in the bound state 
is described by the Schrodinger equation: 



(-Vl + me /h 2 )(f>o = 0, 



(7) 



where m is the atom mass, and the (negative) energy of 
the molecular state is E = —Eq. A general solution of 
Eq. |J7Jl at any finite r can be written as 



Mr) = V'o(r) + /i o G eo (r,0), 



(8) 



where G eo (r, r') = (l/47r|r — r'|) exp(— y meo/h |r — r'|) 
is the Green function of Eq. J7J), and ipo is the solution 
of this equation that is finite and regular at any distance 
including r = 0. For r — > 0, the wavefunction (jHJ should 
satisfy the boundary condition (JTJ- 

One can easily see that any non-trivial solution of 
Eq. J7J, finite at r — > oo, behaves as 1/r for r — > 0. 
Therefore, we have ipo = 0. Then, using the boundary 
condition Q for <f>o(r) at vanishing r, one finds that the 
binding energy of the weakly bound state is 



fe2 / 2 

£o = n /ma , 



(9) 



3 



and the wavefunction </>o normalized to unity is given by 



Mr) 



'2nar) exp(— r/a). 



(10) 



Note that under the condition \a\ 3> R e the main contri- 
bution to the normalization integral comes from distances 
r ^> R e , where Eq. (|10(l is valid. Relative corrections to 
the binding energy e are of the order of (R e /a). 



III. GENERAL FORMALISM FOR THREE 
FERMIONS 

Theoretical studies of the 3-body problem have a long 
prehistory (see |42l | for a review). In this Section we 
consider a three-body system consisting of two identi- 
cal fermions interacting with a third one, which is not 
identical to the first two, via a short-range pair isotropic 
potential U (r) . The fermions have the same mass and we 
will denote the identical ones by the symbol | and the 
third one by the symbol J.. In the center of mass refer- 
ence frame the state of the system with total energy E 
is described by the Schrodinger equation 



-Vt - Vi 



± 



u 



E 



* = Q, (11) 



where y is the distance between the identical fermions, 
and -\/3x/2 is the distance between their center of mass 
and the third atom. Hereinafter we use notations in 
which h — m = 1 . 

The wavefunction \l/ is antisymmetric with respect to 
the permutation of identical fermions, i.e. it changes sign 
under the transformation y — > — y. 

We will discuss the case of resonant two-body interac- 
tion, that is, we assume that the two-body problem for 
the interaction potential U (r) is characterized by a large 
scattering length 



> Re 



(12) 



As has been shown by Efimov 43], in this case short- 
range physics is not important and the 3-body problem 
is universal in the sense that it can be described in terms 
of the 2-body scattering length. One can then use the 
zero range approximation for the interatomic potential. 
This was first done even earlier by Skorniakov and Ter- 
Martirosian 01 i n relation to the neutron-deutron scat- 
tering, which is similar to elastic scatte ring of atoms by 
weakly bound dimers. The work of Ref. [44| was followed 
by related discussions [45j. In this Section, we present 
the form of the zero range approximation for 3-body sys- 
tems, that was outlined in Refs. |33, |46j for the 3-body 
recombination and atom-dimer scattering. 

Under the condition l|12|l . the zero range approxima- 
tion is applicable even at interparticle distances much 
smaller than \a\, as long as these distances greatly exceed 
R e . Also, this approximation properly describes weakly 



bound states of two particles at a > 0. According to the 
zero-range approximation, Eq. I|ll|) is equivalent to the 
Poisson equation 



V 2 + E\* = 0, 



(13) 



with the boundary condition set for a vanishing dis- 
tance between any of the two distinguishable fermions, 
i.e. for r± = (\/3x±y)/2 — * 0. Taking into account the 
symmetry we can write the boundary conditions at the 
two boundaries as 



4tt v ' \r± a 



r± 



0, 



(14) 



where r = =F2y/v / 3 for r± — > 0. The function / contains 
information about the relative motion of a "["-fermion with 
respect to the two other atoms when they are on top 
of each other. In the case of atom-dimer scattering the 
function / plays a role of the wavefunction of the atom- 
dimer relative motion. 

The Green function of Eq. (|13|l . that is the solution 
of this equation with the right hand side equal to <5(x — 
x')<S(y - y'), is given by 



G E {X) 



-{8tt 3 X 2 )- 1 EK 2 (V :Z EX), E < 
i(16w 2 X 2 )- 1 EH 2 {VEX). E>0 



(15) 



Here X = -^/(x — x') 2 + (y — y') 2 , K 2 is an exponen- 
tially decaying Bessel function, and H 2 is a Ha nkel func- 
tion representing an outgoing wave. For -\/|i?|A <C 1 wc 
have 



G = (4^ 3 A 4 )- 1 . 



(16) 



As in the two-body case described by Eqs. I@J and JHJ, 
a general solution of Eq. I|13l) at finite distances between 
t and I fermions can be expressed through the Green 
function (7f H15J) . with coordinates x',y' corresponding 
to a vanishing distance between distinguishable fermions, 
that is for r', — » and for r'_ — > 0. We thus have 



*(x,y) = *o(x,y 
Ge 



/2) 2 + (y + V3r'/2) 2 U(r')dV 



G, 



; (x - r'/2) 2 + (y - V3r'/2) 2 ) h(r')d 3 r', (17) 



where is a properly symmetrized and finite solution 
of Eq. I)13|l . regular at vanishing distances between f and 
I fermions. For a negative total energy E, non-trivial 
solutions of this type do not exist and we have to put 
= 0. The function h(r) has to be determined relying 
on the boundary conditions l|14|) . 

For this purpose we consider the limit r + — > 
and analyze the leading behavior of the terms on the 
right hand side of Eq. (|17fl . The argument of the 
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Green function in the third term can be written as 



rr 



\/3: 



and this term is finite for ra 



_, where r = (x- \/3y)/2, 
-> 0. It can be written as 

(18) 



G E Wr 2 + r' 2 - rr')h(r')d 3 r'. 
In the second term, the argument of the Green function 



takes the form y (r — r') 2 + r^_. We then subtract from 
this term and add to it an auxiliary quantity 



Hr) J G E Lf{ 



d J r 



= (l/47rr + )/i(r)exp(-v /Z £>+). (19) 
The result of the subtraction gives a quantity I = 
J G E (^J(r -r') 2 + r\)[h{Y r ) - h{r)]d 3 r'. For r+ -> 0, 
this quantity remains finite and can be written as 

I = P[ G E (\r - v'\)[h(r') - h{v)]d 3 r\ (20) 

where the symbol P denotes the principal value for the 
integration over dr' . The derivation of Eq. i|20|) and the 
proof of the convergence of the integral on the right hand 
side of this equation are given in Appendix. 

Neglecting terms that are vanishing for r+ — > 0, the 
two last lines of Eq. (|17fl are given by the sum of 
Eqs. JTHJl, (113, and Trm s, in the limit r+ -> 

we can write Eq. 117|l in the form containing a singular 
term proportional to r^ 1 , and regular terms independent 
of r i : 



47T 



-Mr) 



{G E {\v-v'\)[h{v')-h{v)\ 



- G E Wr 2 + r' 2 - rr')/i(r')}d 3 r' + D(t)/4w, (21) 

where D(r) = 47r* (-r/2, \/3r/2), and = -i-y/B 

for positive energies E. Hereafter we omit the principal 
value symbol P for the first term of the integrand of 
Eq. HU). 

As Eq. (J2IJ should coincide with Eq. (|T1J| at r + — > 

0, comparing the singular terms of these equations we 
immediately find that h(r) = /(r). Matching the regular 
terms yields the equation for the function /: 



-£)/(r) = £>(r), 



(22) 



where the integral operator L E is given by 

WW = 47r| {G B (|r-r'|)[/(r)-/(r')] 

+ G^^ + r' 2 -rr')/(r')}dV. (23) 

The operator on the left hand side of Eq. 1221) con- 
serves angular momentum. Therefore, one can expand 
/ and D in spherical harmonics and work only with a 
set of uncoupled equations for functions of a single vari- 
able r. The knowledge of the function /(r) allows one 
to calculate all scattering parameters. In the next sec- 
tion we demonstrate this by calculating the atom-dimer 
scattering length in the ultracold limit. 



IV. ATOM-DIMER ELASTIC SCATTERING 

According to the discussion in Section [HI for a large 
positive scattering length a^> R e one has a weakly bound 
molecular state of two distinguishable fermions (f and |), 
with the binding energy £o = 1/a 2 and the wavefunction 
4>q described by Eq. (|1U|) . We now consider elastic col- 
lisions of these bosonic dimers with f (or J.) fermionic 
atoms at collision energies e < eq. Then the total energy 
of the three-body atom-dimer system is E = (e — £q) < 
and collisions do not lead to dissociation of the dimers. 

From Eq. (|10|l one finds that the size of the weakly 
bound molecular state is ~ a. For a large separation 
between the atom and the dimer, the three-body wave- 
function factorizes into a product: 



*(x, y) « o (r + )x(r), r > a. 



(24) 



A characteristic value of the distance r+ between atoms 
in the dimer is ~ a, and in the limit r 3> a the atom- 
dimer separation is equal to The wavefunction 
x(r) describes the relative atom-dimer motion and can 
be represented as a superposition of an incident and scat- 
tered waves. For r — > oo we have: 



x(r) ps exp(ifcz) + 



2F(k,8) exp(ifer) 
a/3 r ' 



(25) 



where z is the direction of incidence, 8 is the scattering 
angle, the relative wavevector is defined as 2fc/-\/3, and 
F(k, 6) is the scattering amplitude. 
For r + — > we have 



M r +) 



Then, comparing l|24|) with Eq. I|14|) we obtain a rela- 
tion between the functions x an d /: 



X (r) = vV^/W 



(26) 

As we commented in Section IHII for E < we have 
to put = in Eq. |JT7|). This leads to D(r) = 0, and 
Eq. JUJ) takes the form 



-E)f(r) = 0. 



(27) 



The links of our approach to the method of Skorniakov 
and Ter-Martirosian 0] are seen from the fact that 
Eq. H27JI leads to the same equation for the Fourier trans- 
form of the function /(r) as in Ref. [44| . 

We first demonstrate a general approach for solving 
Eq. H27|l and finding the scattering amplitude. For r ^> a, 
the second term on the right hand side of Eq. Q23JI is 
exponentially small and can be omitted. Performing the 
integration in the first term we reduce Eq. (|23|l to the 
form 



L E f{v) = J d 3 q{^-E + q 2 
x /(q)exp(iqr); r > a, 



(28) 
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where /(q) is the Fourier transform of the function /(r). 
Acting twice with the operator (Le + \f—E) on the func- 
tion /(r) and using Eq. (|2*7|) we then obtain 



-E) 2 f{v) 



d 3 q(q 2 - E)f(qyv 



/(r) 



The total energy is given by 

E = -1/a 2 



(29) 



(30) 



and Eq. I|29[l immediately transforms into the equation 
for the free relative atom-dimer motion; 



fc 2 )/(r) = 0; r > a. 



(31) 



The expansion of the function /(r) and the amplitude 
F(k,9) in Legendre polynomials reads (see Appendix): 



/(r) = ]T* J (2Z + l)P/(cosf?)/ ; (r), (32) 

1=0 

oo 

F(k,6)=Y^(2l + l)Pl(coB6)F,(k). (33) 



1=0 



For the function fi(r), which describes the scattering 
with orbital angular momentum I, the superposition of 
an incident and scattered waves satisfying Eq. (|31|l can 
be written as 



Va/8%fl(r) = y/n/2kr{J l+1/2 {kr) 

+ t(2k/V3)F l {k)H l+1/2 {kr)}; r » a, 



(34) 



where H i+1 / 2 is a Hankel function representing an outgo- 
ing wave, and Ji+1/2 is a Bessel function. For kr — > 00 the 
Hankel function is Hi +l / 2 ~ (— y/2/irkr exp(ikr). 
Then, multiplying both sides of Eq. by i l (2l + 

l)p(cos#), making a summation over I, and taking into 
account Eqs. ((23), G3 and Ip^. we arrive at Eq. (J2SJ. 

Using Eq. I|32l) one reduces Eq. (|27|) to a set of uncou- 
pled integral equations for the functions fi(r): 



-E)fi(r)=0; 1 = 0,1,2,. 



(35) 



where the integral operator L l E for a given I is obtained 
by integrating Le over dOr/Air with the weight P;(cos 9). 
Partial scattering amplitudes Fi(k) are then found by 
solving Eqs. I|35|) and fitting the obtained fi(r) with the 
asymptotic expression (|34() . Then, equations l|32|) and 
B33JI give the function /(r) and the total scattering am- 
plitude F(k,9). 

In the ultracold limit, where the condition 



(36) 



is satisfied, the scattering is dominated by the s-wave 
contribution and can be analyzed setting k — and writ- 
ing the total energy as E = —1/a 2 . This is clearly seen 
from Eq. I|34|l with 1 = 0, which for k — > reads: 



/o(r) w v / 87i>(l - 2a Qd /V3r); r > a, 



(37) 



where a Q( i = — -Fo(O) is the atom-dimer scattering length. 

Thus, for finding a a d one should solve Eq. I|35|l with 
I = 0, assuming the limit k — > 0. This equation then 
reduces to 



L° E fo(r) = 0. 



(38) 



Using Eq. I|15|l and integrating over dO r /An in Eq. i|23|) 
we represent Eq. I|38|l in the form: 



L° E fo(r) 



1 



(/o(r)-/o(r')){ 
}+2/o(r'){ 



Kidr'-rl/a) 



irar J 

K x {{r + r')/a)\ , r K^Vr 2 + r' 2 -rr'/a) 



K\ (y/r 2 + r' 2 + rr' / a) 



\/r 2 + r' 2 + rr 



•'/a) \ 
1 J 



r dr 



--0, 



(39) 



where i^i is an exponentially decaying Bessel function. 
It is easily seen that the principal value of the integral in 
the first line of Eq. I|39|l is finite. For r' — > r the integrand 
behaves as l/(r — r'). 

We numerically solved Eq. and found the function 
/o(r) at all distances r. Fitting the obtained fo{r) at 
r a with the asymptotic expression l|37l) we arrived 
at the atom-dimer scattering length a a d = 1.2a, which 
reproduces the result of Ref. 01 • 0m calculations also 
show that the behavior of fo suggests a soft-core atom- 
dimer repulsion, with a range of the order of a. 



V. RELAXATION IN ATOM-DIMER 
COLLISIONS 

Weakly bound dimers that we are considering are 
diatomic molecules in the highest ro-vibrational state. 
Hence, they can undergo relaxation into deep bound 
states in their collisions with each other or with unbound 
atoms. The released binding energy of a deep state is 
~ h 2 /mR 2 . It is transformed into the kinetic energy of 
particles in the outgoing collisional channel and they es- 
cape from the trapped sample. Therefore, the process of 
collisional relaxation of weakly bound dimers determines 
the lifetime of a gas of these molecules and, in particular, 
possibilities to Bose-condense such a gas. 

In our previous work |3S| we have shown that this pro- 
cess is suppressed due to Fermi statistics for the atoms 
in combination with a large size of the dimer. The bind- 
ing energy of the dimers is £0 = Ti 2 /ma 2 and their size 
is ~ a 3> R e - The size of deep bound states is of the 
order of R e . Hence, the relaxation requires the presence 
of at least three fermions at distances ~ R e from each 
other. As two of them are necessarily identical, due to 
the Pauli exclusion principle the relaxation probability 
acquires a small factor proportional to a power of (qR e ), 
where q ~ 1/a is a characteristic momentum of the atoms 
in the weakly bound molecular state. 

In this Section we discuss relaxation of weakly bound 
dimers into deep bound states in ultracold atom-dimer 



6 



collisions satisfying the condition i|36|) . Relying on the 
inequality a ^> R e we develop a method that allows us 
to establish a dependence of the relaxation rate on the 
scattering length a, without going into a detailed anal- 
ysis of the short-range behavior of the systems of three 
atoms. It is assumed that the amplitude of the inelastic 
process of relaxation is much smaller than the amplitude 
of elastic scattering. Then the dependence of the relax- 
ation rate on a is related only to the a-dependence of the 
initial-state 3-body wavefunction '5. 

The relaxation occurs when all of the three atoms ap- 
proach each other to distances of the order of R e . At 
these interatomic distances, as well as at all distances 
x, y <C a, the wavefunction W in the ultracold limit is de- 
termined by the Schrodinger equation with E = 0. 
Therefore, it depends on the scattering length a only 
through a normalization coefficient: 



* = A(a)ip; x,y<a, 



(40) 



where the function ip is independent of a. The probability 
of relaxation and, hence, the relaxation rate constant a a d 
are proportional to 1^1 2 at distances x,y ~ R e . We thus 
have 



a ad cx |A(a)| 5 



(41) 



The goal then is to find the coefficient A(a), which de- 
termines the dependence of the relaxation rate on a. 

In the region where R e <C {x, y} <C a, the a- 
independent function i\) can be found in the zero-range 
approximation. Then, matching the wavefunction \E' 
given by Eq. I|4(J|) . with the result of the zero range ap- 
proximation at interparticle distances larger than a gives 
the coefficient A(a). 

We start with analyzing the behavior of the 3-body 
wavefunction in the zero range approximation at dis- 
tances x,y -C a. In this region of distances, $ is recon- 
structed through the function /(r) from Eq. I|17|) with 
h(r) = /(r) and "to = 0, setting E — > 0. Accordingly, 
one should use the Green function Go <|16[l m this equa- 
tion, which then reads: 



*(x,y) 



r'/2) 2 + (y + V3r'/2y 



Go(v / (x-r'/2)2 + ( y -v / 3r'/2)2) 
x /(r')dV; x, y < a. 



(42) 



The main contribution to the integral in Eq. I|42|l comes 
from distances r' <C a and, hence, we have to find the 
function / at these distances. In the ultracold limit where 
the condition (|36|) is satisfied, the inelastic process of 
relaxation is dominated by the contribution of the s-wave 
channel. As we found in Section llVl the s-wave part /o of 
the function / is determined by Eqs. I|38|l and (|39|l . For 
r <C a, the distances r' in Eq. (|39|l are also much smaller 
than a. This is equivalent to setting the limit E — > 



(a 1 — > 0) for the integral operator L° E in Eq. (|39|l . and 
this operator then reduces to 



L° fo(r) 



1 



{(r' - 



(/o(r) - /o(0) 

■2 



(r' — r) 2 (r + r') 2 
1 



l - } 

12 i rr l J 



(43) 



where the integration of the term containing (r — r') 2 
assumes a principal value of the integral. 

Thus, the function fo(r) at r < o is a solution of the 
integral equation 



L° fo(r) = 0. 



(44) 



The operator Lq has a property that L^r" = A(^)r ly_1 
for —5 < Re (v) < 3, and the integral in Eq. (|43|l diverges 
outside this interval. The function \{y) is given by 



\(y) = (v + l)tan— + 



ttv 4 sin[7r(;y + l)/6] 



V3 cos(?r^/2) 



(45) 



In the specified interval of v the function A has two roots, 
v + = 1.166 and v_ = —3.166. Accordingly, the solution 
of Eq. 14411 is a linear superposition: 



fo(r) « C + r v + + C_r v 



r <C a. 



(46) 



The determination of the ratio G + /G_ involves short- 
range physics and is beyond the scope of this paper. 
However, in the absence of a three-body resonance the 
matching procedure implies that at distances r ~ R e 
both terms in Eq. (|46|l are of the same order of magni- 
tude. So C+/C- cx R~ 4 332 , and at distances r ^> R e one 
has fo(r) w C+r^ . Substituting this result into Eq. l(l2")l 
we find that at distances x,y <C a the three-body wave- 
function takes the form 

*«$°(r!)/(p)/p = G + a> (r!)^- 1 ; x,y<a, (47) 



where p = \/ x 2 + y 2 is the hyperradius, and the set of hy- 
perangles denotes all the other coordinates. Although 
one can explicitly write down the function $°(f2), for us 
it is only important that this function is a-independcnt. 

Comparing Eq. I|47l) with Eq. (|40l) we see that one may 
set ip = &°(£l)p"+~ 1 in the interval of distances where 
R e <C {x,y} <C a. We then have C + = A(a), i.e. the 
function /o can be written as 



fo(r) w A( a y 



r <C a. 



(48) 



Numerical integration of Eq. I|39|) shows that /o smoothly 
interpolates between the asymptotic expressions given by 
Eq. |gU for r < a, and by Eq. (J37JI for r > a. This 
procedure provides matching of the two asymptotes at 
r ~ a and gives the coefficient 



A(a) 



oc a 



(49) 
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In fact, the result of Eq. I|49(l can be obtained in a 
more elegant way, using only Eqs. I|37|) and In the 

zero range approximation the only distance scale is the 
2-body scattering length a. Hence, we may rescale the 
coordinate and represent the function /o in the form 



/o(r) = B(a)f (r/a), 



(50) 



where /o depends on a only through the rescaled coordi- 
nate r/a. The coefficient B(a) is independent of the co- 
ordinate and can be obtained by comparing Eq. I|50|> with 
asymptotic expressions (|3"7l and (|I5|l . Using Eq. i(3"Y|) we 
see that B oc a -1 / 2 , whereas the comparison of Eq. I|50(l 
with Eq. igBJl gives B oc A(a)a l ' + . This immediately leads 
to Eq. P^l) for the coefficient A(a). 

As the dependence of the relaxation rate constant a ad 
on the 2-body scattering length is governed by Eq. 1)41(1 . 
using Eq. we obtain a ad oc a~ s , where s = 1 + 
2v + ~ 3.33. The absolute value of the relaxation rate is 
determined by the contribution of interparticle distances 
~ R e , where the zero-range approximation is not valid. 
This approximation only gives a correct dependence of 
the relaxation rate on a. 

Assuming that the short-range physics is characterized 
by the length scale R e and by the energy scale ft 2 /mR 2 , 
we can restore the dimensions and write the following 
expression for the rate constant of relaxation of weakly 
bound dimers into deep bound states in atom-dimer ul- 
tracold collisions: 



ld = C(hR e /m)(R e /aY 



3.33. 



(51) 



One clearly sees that the relaxation rate rapidly decreases 
with increasing the 2-body scattering length. However, 
the coefficient C depends on a particular system. 

The relaxation due to atom-dimer scattering with 
non-zero orbital angular momenta is very small. Since 
the atom-dimer effective interaction has a characteristic 
range ~ a (see Sec. II VL the p-wave part of the 3-body 
wavefunction "F at short interparticle distances is pro- 
portional to ka. Hence, the p-wave contribution to the 
relaxation rate is oc (ka) 2 and can be omitted for ultra- 
cold collisions satisfying the condition (|3"rJ)l . 



VI. ELASTIC DIMER-DIMER SCATTERING 

As we already mentioned in Introduction, elastic in- 
teraction between weakly bound bosonic molecules of 
j and I fermionic atoms is important for understand- 
ing the physics of their Bose-Einstein condensation and 
for studying the BCS-BEC cross-over in two-component 
atomic Fermi gases. The dimer-dimer elastic scattering is 
a 4-body problem described by the Schrodinger equation 



{ 



V 2 



R 



+ U(n) + U{r 2 ) 



^C/[(r 1 +r 2 ±V2R)/2]-^}* = 0. (52) 



Here we again use units in which % = m = 1. The 
distance between two given f and j fcrmions is r^, and 
r 2 is the distance between the other two. The distance 
between the centers of mass of these pairs is R/^/2, and 
(ri + r 2 ± \/2R)/2 are the separations between ] and [ 
fcrmions in the other two possible || pairs. The total 
energy is E = — 2eq + £, with e being the collision energy, 
and £o = — 1/a 2 the binding energy of a dimer. 

The wavefunction 'F is symmetric with respect to the 
permutation of bosonic jj. pairs and antisymmetric with 
respect to permutations of identical fermions: 

*(ri,r 2 ,R) = <F(r 2 ,n,-R) 

/ 1-! + r 2 ± V2R ri+rgJF %/2R ri -r 2 \ 
{ 2 ' 2 ' V2 J' 

(53) 

The weak binding of atoms in the dimer assumes that 
the 2-body (positive) scattering length is a ^> R e , and 
we employ the zero range approximation in our analysis 
of the molecule-molecule scattering. This is done relying 
on the formulation of this approximation given in Sec. lIIII 
for the three-body problem. Thus, the 4-body system is 
described by the free-particle Schrodinger equation 



[V 2 i+ V 2 2 +VR + i?]* = 0, 



(54) 



and the 4-body wavefunction ^> should satisfy the Bethe- 
Peierls boundary condition for a vanishing distance in 
any pair of "f and J, fcrmions, i.e. for ri — *■ 0, r 2 — > 0, and 
ri + r 2 ± \/2R — > 0. Due to the symmetry it is necessary 
to require a proper behavior of ^ only at one of these 
boundaries. For ri — > the boundary condition reads: 



*(n,r 2 ,R) -> /(r 2 ,R)(l/47iTi - 1/4^) 



(55) 



The function /(r 2 ,R) is analogous to that defined in 
Sec. IIIII and it contains the information about the sec- 
ond pair of particles when the first two are sitting on top 
of each other. 

In the ultracold limit, where the condition l|3()|) is sat- 
isfied, the scattering is dominated by the contribution of 
the s-wave channel. The inequality l(36|l is equivalent to 
£ -C £o and, hence, the s-wave scattering can be analyzed 
from the solution of Eq. ifoU) with E = — 2e < 0. For 
large R the corresponding wavefunction is given by 

* « (fo(n)to(ra)(l - V2a dd /R); R > a, (56) 

where a d d is the dimer-dimer scattering length, and the 
wavefunction of the weakly bound dimer is given by 
Eq. HlOjl. Combining Eqs. (|53|l and (|5*6*)l we obtain the 
asymptotic expression for / at large distances R: 

/(r 2 , R) pa (2/r 2 a) exp (-r 2 /a)(l - V2a dd /R); R > a. 

In the case of the s-wave scattering the function / de- 
pends only on three variables: the absolute values of r 2 
and R, and the angle between them. We now derive and 
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solve the equation for /. The value of add is then deduced 
from the behavior of / at large R governed by Eq. I|57l) . 

We first establish a general form of the wavefunction 
iff satisfying Eq. (|54|l . with the boundary condition l|55|) 
and symmetry relations (|53|l . In our case the total energy 
E = — 2/ a 2 < 0, and the Green function of Eq. (f54"j> reads 

G{X) = (27r)- 9 / 2 (Xa/V2)- 7 ^K r/2 (V2X/a), (58) 

where X = \S - S'\, and S = {ri,r 2 ,R} is a 9- 
component vector. Accordingly, one has \S — S'\ — 
v/(ri - r'i) 2 + (r 2 - r' 2 ) 2 + (R - R') 2 - In analogy with 
the 3-body case, the 4-body wavefunction iff can be ex- 
pressed through G(\S — S'\) with coordinates S' cor- 
responding to a vanishing distance between f and J, 
fermions, i.e. for r'i — > 0, r' 2 — > 0, and (r'i + r' 2 ± 
\/2R')/2 — > 0. Thus, for the wavefunction satisfying the 
symmetry relations 1)530 we have 

= *o+ / dVd 3 i?'[G(|S-Si|) + G(|S , -S 2 |) 



(59) 



-G(|5-5+|) -GQS-S-I) /i(r',R) 



where Si = {0,r',R'}, S 2 = {r',0,-R}, and 5± = 
{r / /2±R'/v / 2, r'/2TR7v / 2, Tr'V2}. The function * is 
a properly symmetrized finite solution of Eq. 1)541) , regular 
at vanishing distances between | and j fermions. For 
i? < 0, non-trivial solutions of this type do not exist and 
we have to put iffo = 0. The function h(r 2 ,R) has to 
be determined by comparing iff (|59|l at ri — > 0, with the 
boundary condition (J5SJ. 

This procedure is similar to that developed in Sec- 
tion lllll for the 3-body case. Considering the limit ri — > 
we extract the leading terms on the right hand side of 
Eq. (|59l) . These are the terms that behave as l/r\ or 
remain finite in this limit. The last three terms in the 
square brackets in Eq. I|59|) provide a finite contribution 



d 3 r'd 3 R'h(r',R') G{\S 2 - S 2 \) 



- G(\S 3 - S+\) - G(\S 2 - S-\) 



(60) 



where §2 = {0,r2,R}. For finding the contribution of 
the first term in the square brackets we subtract from 
this term and add to it an auxiliary quantity 



h(r 2 ,R) J GQS-S 1 \)dPr f dPB! = 



h{r 2 ,R ) 



■ exp ( 



-V2~n/a). 



(61) 



The result of the subtraction yields a finite contribution 
which for n — » can be written as 

J d 3 r'd 3 R'[h(r',R') - h(r 2 ,R)]G{\S - Si\) = 



P / d 3 r , d 3 R'[h(r' ,R') - h(r 2 ,R)]G(\S 2 - Si\ 



with the symbol P standing for the principal value of the 
integral over dr' (or dR'). Equation (|6^)l is derived in 
Appendix and it is proven that the integral in the second 
line of this equation is convergent. 

In the limit r% — * 0, the right hand side of Eq. I|61|l is 
equal to 



h(r 2 ,R)(l/47rri - y/2/Aita). 



(63) 



We thus find that for ri 
Eq. lt59"|) takes the form 



the wavefunction iff of 



*(ri,r 2) R) 



h{r 2 ,R) 



'fnr-i 



K; ri->0, (64) 



where 1Z is the sum of regular n -independent terms given 
by Eqs. I|60() and (|62H . and by the second term on the 
right hand side of Eq. l|63Jl . Equation Ij64() should coincide 
with Eq. I|55|) , and comparing the singular terms of these 
equations we find h(r 2 , R) = f(r 2 , R). As the quantity TZ 
should coincide with the regular term of Eq. I|55() , equal 
to — /(r 2 , R)/4-7ra, we obtain the following equation for 
the function /: 

d 3 r'd 3 R'{G(\S ~ S^tfiv^R') - f(v,R)} 
+ [G(\S- S 2 |) - S ± |)]/(r',R')} 
= (V2-l)/(r,R)/47ro. (65) 



Here S — {0, r, R}, and we omitted the symbol of prin- 
cipal value for the integral in the first line of Eq. H65|) . 

As we already mentioned above, for the s-wave scat- 
tering the function /(r, R) depends only on the absolute 
values of r and R and on the angle between them. Thus, 
Eq. I|65|) is an integral equation for the function of three 
variables. We have solved this equation numerically for 
all distances R and r, and all angles between the vec- 
tors R and r. Fitting the asymptotic expression (|57ll at 
R 3> a with the function /(r,R) obtained numerically 
from Eq. H65|) , we find with 2% accuracy that the dimer- 
dimer scattering length is: 



add = 0.6a > 0. 



(66) 



(62) 



Our calculations show the absence of 4-body weakly 
bound states, and the behavior of / at small R suggests 
a soft-core repulsion between dimers, with a range ~ a. 

The result of Eq. 166|) is exact, and its derivation was 
outlined in our previous work [38j. Equation l|66|) in- 
dicates the stability of molecular BEC with respect to 
collapse. Compared to earlier studies which assumed 
add = 2a 4], Eq. (|66|l gives almost twice as small a 
sound velocity of the molecular condensate and a rate 
of elastic collisions smaller by an order of magnitude. 
We should mention here that the result of earlier studies 
was reconsidered in Ref. by using a diagram- 
matic approach which leads to add = 0.75a. However, 
this approach misses a number of diagrams which give 
a contribution of the same order of magnitude as those 
taken into account. 
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VII. RELAXATION IN DIMER-DIMER 
COLLISIONS 

In this Section we generalize the results obtained in 
Sec. to the relaxation of weakly bound dimers into 
deep bound states in dimer-dimer collisions. We again 
consider the ultracold limit described by the condition 
(|36[) . where the relaxation process is dominated by the 
contribution of the s-wave dimer-dimer scattering. The 
key point of our discussion is that the dimer-dimer re- 
laxation collisions are to a large extent similar to the 
atom-dimer ones. 

Indeed, the relaxation process requires only three 
atoms to approach each other to short distances of the or- 
der of R e . The fourth particle can be at a large distance 
from these three and, in this respect, does not partici- 
pate in the relaxation process. This distance is of the 
order of the size of a dimer, which is ~ a ^> R e . As well 
as in the case of atom-dimer collisions, the dependence 
of the relaxation rate on the 2-body scattering length 
a is determined by the a-dependence of the initial state 
wavefunction "P. We thus see that the configuration space 
contributing to the relaxation probability can be viewed 
as a system of three atoms at short distances ~ R e from 
each other and a fourth atom separated from this sys- 
tem by a large distance ~ a. In this case the 4-body 
wavefunction decomposes into a product: 



(67) 



where ^P^ 3 ^ is the wavefunction of the 3-fermion system, 
p and f2 are the hyperradius and the set of hyperangles 
for these fermions, z is the distance between their cen- 
ter of mass and the fourth atom, and the function 7^(z) 
describes the motion of this atom. Note that Eq. 1)6 7|) 
remains valid for any hyperradius p <C |z| ~ a. 

The transition to a deep bound 2-body state occurs in 
the system of three atoms and does not change the wave- 
function of the fourth atom, rj(z). Therefore, averaging 
the transition probability over the motion of the fourth 
particle, the rate constant of relaxation in dimer-dimer 
collisions can be written as 



a ( id 



,(3) 



\T 1 {z)\ 2 d- i z = a 



(3) 



(68) 



where is the rate constant of relaxation for the 3- 
atom system. 

We thus obtain that the problem is reduced to the 
relaxation in atom-dimer collisions. The difference from 
the case discussed in Section is that now the relative 
momentum of the collision is ~ 1/a. This is seen by 
considering at large distances between the dimer 

and the fermionic atom, hereafter referred to as the 3-rd 
fcrmion. This fermion is in the bound molecular state 
with the fourth atom. As the size of this state is ~ a, 
the expansion of the wavefunction of the 3-rd fermion in 
plane waves shows that its momentum q is of the order 
of 1/a. 



Keeping in mind the discussion in Section[V]we see that 
the result for the relaxation rate following from Eq. (|5 1|) 
remains valid for the dimer-dimer collisions. The fact 
that the relative momentum is q ~ 1/a can only change 
the numerical coefficient, not the dependence of the re- 
laxation rate constant on the 2-body scattering length a. 
However, the result of Eq. I|51|l accounts only for the s- 
wave scattering of the 3-rd fermionic atom on the dimer, 
which provides the leading relaxation channel for ultra- 
cold atom-dimer collisions. In the case of dimer-dimer 
collisions, there is a relaxation channel that is more im- 
portant in the limit of large a (see below). For the s- 
wavc dimer-dimer scattering, both the 4-th and the 3-rd 
fermions (bound to each other in the molecular state) can 
undergo the p-wave scattering on the other dimer in such 
a way that their total orbital angular momentum is equal 
to zero. We thus should consider the relaxation for the 
p-wave collisions between the 3-rd fermionic atom and 
the dimer. These collisions are not suppressed as their 
relative momentum is q ~ 1/a. 

As we discussed in Section |VJ for the hyperradius 
p -C a the wavefunction of the 3-fermion system, ty( 3 ', 
depends on the scattering length a only through a nor- 
malization coefficient and can be written in the form 14011 . 
We thus have Vl^ 3 ) = A(a)i/j, with a-independent func- 
tion ip 7 and the relaxation rate constant depends on a 
as add oc |^4(a)| 2 . This requires us to find the coefficient 
A(a) for the case of the p-wave scattering of the 3-rd 
fermionic atom on the molecule. We will do this in the 
zero range approximation, along the lines of our discus- 
sion of the s-wave atom-dimer scattering in Section Ivl 

We first consider the region of interparticle distances 
where R e <C p <C a. Then the function /(r) for the 
3-body problem is determined by Eq. (|27|) in which we 
have to set the limit E — > and a -1 — ► 0. Therefore, 
the integral equation l|35|l for the p-component of the 
function / takes the form: 



(69) 



where the operator Lq is obtained substituting the ex- 
pansion (|32|l into Eq. i|23|) with the Green function Go 
(|16|l . multiplying by cos 9, and integrating over the an- 
gles. This yields 



1 f°° 
Llhir) = - I dr 



'(7[ (/l(r) - /l(/)) {f^ 

(r + r') 2 } + 2/l(r,) {^T7 



r — r) 
1 



1 



(r + r') 



M 2 ) 



2/i(r') l ^{r + r'){r 2 + r' 2 - rr') 
y/\r — r'|(r 2 + r' 2 + rr') 



(70) 



where one should take a principal value for the integral 
of the term containing (r — r 1 ) 2 . 

Like the operator Lq, the operator Lj has a property 
L\r v = X^r"- 1 for -4 < Re(i/) < 2. The function 
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\{y) is now given by 

4 vcos [k(v + l)/6] — 2sin(7ri//6) 



A(") = 



(z/ + 1) sin(7ri//2) 



fff + 2) 7W 

cot — 

i/+ 1 2 



(71) 



and has two roots: v + = 0.773 and v_ = —2.773. Then, 
with the same arguments as in between equations (|45f) 
and H48fl in Section [V] we obtain short-distance expres- 
sions for the function /i and for the part of the 3-body 
wavefunction corresponding to the p-wave scattering 
of the 3-rd fermion on the molecule: 



f 1 {r)ttA(a)a v +{r/aY+, 
tfPO = A(a)$ l {Q)p u +- 1 , 



r <C a; 
p <C a, 



(72) 
(73) 



where the function is independent of a. 

For relative momenta k ~ 1/a, the 2-body scattering 
length remains the only distance scale of the zero range 
approximation. In particular, the p-wave scattering am- 
plitude in Eq. (|34|l will be F\ ~ a. We therefore can write 
the function /i in the form lj5U|) : /i(r) = B(a)fx(r/a), 
where /i depends on a only through the rescaled coor- 
dinate r/a. In order to be consistent with Eq. (|34|l for 
k ~ 1/a, we have to put -B(a) oc a -1 / 2 . Then the com- 
parison of the resulting /i with Eq. I|72|l gives the coeffi- 
cient A(a) oc a~ 1 / 2 ~ v + and the relaxation rate constant 
is add oc |A|(a)| 2 oc a _s , with ,s = l + 2v + ~ 2.55. Restor- 
ing the dimensions we have 



add = C(hR e /m)(R e /aY 



s = 2.55. 



(74) 



One can think of the relaxation mechanism, where the 
scattering of the 3-rd fermion on the dimer occurs with 
higher orbital angular momenta I. In this case the 4-th 
atom scatters on the dimer with the same I, and the total 
angular momentum of the dimer-dimer scattering should 
be equal to zero. Our analysis shows that these scattering 
mechanisms lead to a power law dependence add oc a~ s , 
with larger values of s than in Eq. I|74l) . Hence, for large 
a they can be neglected. 

Equation (|74|l shows a slower decrease of the relaxation 
rate with increasing a than in the case of atom-dimer 
collisions. Obviously, in the limit R e /a — > the dimer- 
dimer relaxation should dominate over the atom-dimer 
one. The competition between these two relaxation pro- 
cesses can be present if the ratio R e /a is not too small 
and the densities of dimers and atoms in the gas are com- 
parable with each other. 



VIII. WIDE AND NARROW FESHBACH 
RESONANCES 

In experiments with alkali atom gases large values of 
the 2-body scattering length are achieved by using Fesh- 
bach resonances. In the vicinity of the resonance, the 
2-body problem is characterized by a strong coupling 



between the continuum states of colliding atoms and a 
bound molecular state of another hyperfine domain of 
these atoms. The resulting scattering length depends on 
the detuning from the resonance, i.e. on the energy dif- 
ference A between the border of the continuum of collid- 
ing atoms and the bound molecular state. The splitting 
between the two hyperfine domains and, hence, the de- 
tuning A depend on the magnetic field, which makes the 
scattering length tunable by varying the field. 

One thus has a two-channel problem which can be de- 
scribed in terms of Breit-Wigner scattering |4ll |48| , the 
open channel being the states of colliding atoms and the 
closed channel the bound molecular state of the other hy- 
perfine domain. Various aspects of this type of problems 
have been discussed by Feshbach and Fano [H(J. In 
cold atom physics the idea of Feshbach resonances was 
introduced in Ref . jHJ , and optically induced resonances 
have been discussed in Refs. |52ll53|. 

We now analyze to which extent our results for 3-atom 
and 4-atom systems of fermions with a positive 2-body 
scattering length a ^> R e describe the situation of Fes- 
hbach resonances. For low collision energies e, omitting 
the (small) background scattering length, the scattering 
amplitude is given by |4l| : 



F(s) 



(75) 



where the quantity h^f / ^/m = W characterizes the cou- 
pling between the two hyperfine domains. In Eq. (|75|l the 
detuning A is positive if the bound molecular state is be- 
low the continuum of colliding atoms. Then for A > 
we have a positive scattering length near the resonance, 
a = —F(0) = W/A. Introducing a characteristic length 



R* = h 2 /mW 



(76) 



and expressing the scattering amplitude through the rel- 
ative momentum of particles k — ^me/h, equation Ij75(l 
takes the form: 



F(k) 



1 



+ R*k 2 + ik' 



(77) 



The validity of Eq. I|77|) does not require the inequality 
kR* <C 1. At the same time, Eq. I|77|l formally coincides 
with the amplitude of scattering of slow particles by a po- 
tential with the same scattering length a and an effective 
range R = — 2R* , valid under the condition kR <C 1. 

The length R* is an intrinsic parameter of the Feshbach 
resonance problem. It characterizes the width of the res- 
onance. From Eqs. (|75|l and (|76|l we see that large W and, 
hence, small R* correspond to a wide resonance, whereas 
small W and large R* lead to a narrow resonance. The 
issue of wide and narrow resonances is now actively be- 
ing discussed in literature [3 |H |H |H 111 We 
would like to point out here that the use of terms "wide" 
and "narrow" depends on the problem under considera- 
tion. For example, in the unitarity limit where a — > ±oo, 
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Eq. (|77jl shows that the length R* drops out of the prob- 
lem under the condition kR* C 1. In a quantum de- 
generate Fermi gas the characteristic momentum of par- 
ticles is the Fermi momentum kp — (i^n) 1 ^ , where 
n is the gas density. Thus, the inequality kpR* <C 1, 
referred to as the condition of a wide resonance, en- 
sures universality of the problem 0, 0, [5(| E3, l58| - 
The only length and energy scales in the gas will be the 
mean interparticle separation rT 1 ! 3 and the Fermi energy 
Ep = (3n 2 n) 2 / 3 h 2 /2m, and the system acquires univer- 
sal thermodynamics [lq . 

In our case the situation is different. We are consider- 
ing weakly bound diatomic molecules in the open channel 
with a binding energy @, and our discussion of atom- 
dimer and dimer-dimer collisions assumes that there is a 
weakly interacting gas of these dimers and atoms. The 
most important limitation is related to the binding en- 
ergy and the wavefunction of the dimers. The energy of 
the weakly bound molecular state is determined by the 
pole of the scattering amplitude l|77l) . One then finds 
that this state exists only for a > and the universal 
expression (0 used in our calculations, £o = h 2 /ma 2 , 
requires the inequality [5^| 



R* < a. 



(78) 



Under this condition the wavefunction of the weakly 
bound molecular state has only a small admixture of the 
closed channel |55j . 

Our calculations for the atom-dimer and dimer-dimer 
collisions were done in the ultracold limit where ka <C 
1. In a thermal gas of atoms and dimers the char- 
acteristic momentum is the thermal momentum kx — 
(2mT /Ti 2 ) 1 / 2 , whereas for a degenerate, in particular 
Bose-condensed gas of the dimers, the characteristic mo- 
mentum is the inverse healing length, (na) 1 / 2 . In both 
cases the inequality ka <C 1 assumes the weakly inter- 
acting regime, where na 3 <C 1. One can also see that in 
the limit k«l the inequality Eq. I|78|) makes the scat- 
tering amplitude F(k) (|77() momentum independent and 
equal to —a. This justifies the use of our single-channel 
zero range approximation for calculating atom-dimer and 
dimer-dimer interactions and collisional properties. We 
thus obtain that for our problem the condition of a wide 
Feshbach resonance is given by Eq. (|78|) . Under this con- 
dition the problem is universal in the sense that the size of 
weakly bound dimers, and atom-dimer and dimer-dimer 
scattering properties are characterized by a single param- 
eter, the 2-body scattering length a. 

Note that the interaction between the two channels 
of the Feshbach problem is efficient at interparticle dis- 
tances which are of the order of R e . Therefore, the Fesh- 
bach character of scattering does not influence the condi- 
tion a 3> R e that allows us to use the zero range approx- 
imation for the entire region of interparticle distances. 



IX. CONCLUDING REMARKS 

In most experiments with weakly bound diatomic 
molecules produced by using Feshbach resonances in a 
degenerate two-component atomic Fermi gas, the wide 
resonance condition (jTSj) was satisfied. This was the case 
with 40 K 2 molecules at JILA H3|Ml3l, and with 6 Li 2 
at Innsbruck pHIMIIIOSI, MIT Duke 0H3, 

ENS HHH UK, and Rice [12 In these experimental 
studies the length R* is of the order of or smaller than 
20A, and for the achieved values of the scattering length 
a from 500 to 2000A the ratio R*/a is smaller than 0.1. 
The only exception is the experiment at Rice with 6 Li 
near a narrow Feshbach resonance at 543 G [26j . For this 
resonance the length R* is large and at obtained values of 
a the condition (|78|l is not fulfilled. Therefore, the Rice 
experiment |26| cannot be described by our theory 

Experimental studies of dimers produced in a Fermi 
gas by using wide Feshbach resonances at JILA, Inns- 
bruck, MIT, and ENS are well described within our the- 
oretical approach. It should be mentioned that for a > 
and equal concentrations of the two atomic components 
of the gas, at temperatures well below Ep practically all 
atoms should be converted into dimers if the gas density 
satisfies the inequality na 3 <C 1 59]. In ongoing exper- 
iments the imbalance between the atomic components 
is fairly small, and at sufficiently low temperatures there 
can only be a small fraction of unpaired fermionic atoms. 

The results at JILA_|27 
MIT [M H, ENS E 



33J, Innsbruck [25j,|29, 36], 
31| . and Rice [32 show a re- 
markable collisional stability of weakly bound diatomic 
molecules 40 K 2 and 6 Li2. At molecular densities n ~ 10 13 
cm~ 3 the lifetime of the gas ranges from tens of millisec- 
onds to tens of seconds, depending on the value of the 
scattering length a. A strong decrease of the relaxation 
rate with increasing a, following from Eq. (|74|l . is consis- 
tent with experimental data. The potassium exp eriment 
at JILA [53 and lithium experiment at ENS [31( give the 
relaxation rate constant add oc a~ s , where s w 2.3 with 
15% of accuracy for K 2 , and s ~ 2.0 with 40% of accu- 
racy for Li 2 . The absolute value of the rate constant for 
a thermal gas of Li 2 is add ~ 2 x 10 13 cm 3 /s for the scat- 
tering length a ~ 800 A 24]. For K 2 it is by an order of 
magnitude higher at the same value of a [27]], which can 
be a consequence of a larger value of the characteristic 
radius of interaction R e . 

At realistic temperatures the relaxation rate constant 
add is much smaller than the rate constant of elastic col- 
lisions &ita 2 dd VT, where vt is the thermal velocity. For 
example, for the Li 2 weakly bound dimers at a tempera- 
ture T ~ 3/xK and a ~ 800A, the corresponding ratio is 
of the order of 10" 4 or 10~ 5 . This opens wide possibili- 
ties for reaching BEC of the dimers and cooling the Bose- 
condensed gas to temperatures of the order of its chemical 
potential. Long-lived BEC of weakly bound dimers has 
been recently observed for 40 K 2 at JILA [28Ll33T| and for 
6 Li 2 at Innsbruck [MEl, MIT [MHI, ENSlp, and 
Rice [32 ■ Measurements of the dimer-dimer scattering 
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length in these experiments confirm our result add = 0.6a 
|38| with accuracy up to 30% [6(j ■ This result is also 
confirmed by recent Monte Carlo calculations of the 
ground state energy in the molecular BEC regime. 

In conclusion, we have developed a theory of elastic and 
inelastic collisions of weakly bound molecules formed in 
a two-component atomic Fermi gas. We emphasize that 
the remarkable collisional stability of these molecules at 
a 3> R e is due to Fermi statistics. This effect is not 
present for weakly bound molecules of bosonic atoms, 
even if they have the same large size. Indeed, identical 
fcrmionic atoms participating in the relaxation process 
at short interparticle distances have very small relative 
momenta k ~ 1/a and, hence, the process is suppressed 
compared to the case of molecules of bosonic atoms. 

The long lifetime of weakly bound diatomic molecules 
of fermionic atoms allows interesting manipulations with 
these dimers. It seems realistic to arrange a deep evap- 
orative cooling of their Bose-condensed gas to temper- 
atures of the order of the chemical potential. Then, 
converting the molecular BEC into fermionic atoms by 
adiabatically changing the scattering length to negative 
values, one provides an additional cooling. The obtained 
atomic Fermi gas will have extremely low temperatures 
T ~ 10~ 2 Ep and can be already in the BCS regime 
[62l| . At these temperatures one has a very strong Pauli 
blocking of elastic collisions and expects the collisionless 
regime for the thermal cloud, which is promising for iden- 
tifying the BCS-paired state through the observation of 
collective oscillations or free expansion [63|. 

Another idea is related to transferring weakly bound 
molecules of fermionic atoms to their ground ro- 
vibrational state b y u sing two-photon spectroscopy, as 
proposed in Ref. |6J] for molecules of bosonic atoms. 
Long lifetime of weakly bound dimers of fermionic atoms 
at densities ~ 10 13 cm~ 3 should provide a much more 
efficient production of ground state molecules compared 
to the case of dimers of bosonic atoms, where one has 
severe limitations on achievable densities and lifetimes. 
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Appendix 

We first prove that in the limit r + — > we have Eq. i|20|) 
and the integral on the right hand side of this equation 
is convergent. So, we consider the integral 



G, 



;(^/(r - r') 2 + r\ )[h{v') - h(v)]d 3 r' (79) 



and analyze the limiting case of r + — > 0. For this purpose 
we use the expansion of the functions h(r') and h(r) in 
spherical harmonics. This type of expansion can be made 
for any function of the components of a three-dimensional 
vector and it reads 



h ( r ) = J2 hi m (r)Y lm (9 r ,<p r ), 



(80) 



--1 



with 9 r and <p r being the polar and azimuthal angles of 
the vector r with respect to a quantization axis. Note 
that in the analysis of the atom-dimer elastic scattering 
in Section II VI we used an expansion in Legendre poly- 
nomials for the function f(r). The reason is that the 
bound molecular state has zero orbital angular momen- 
tum and there is a cylindrical symmetry in the system. 
For a general 3-body problem one should, in principle, 
use Eq. (JH01). 

For spherical harmonics Yj m (0 r i , 4> r > ) in the expansion 
(|80|l for the function h(r'), we use the relation (see |4lj): 



Y lm (e r ,,^) =J2D™> m (<t>r,0r,0)Y lm ,(e', ((>'), (81) 



where D^, m is the matrix of finite rotations, and 
6',4>' are the polar and azimuthal angles of the vec- 
tor r' with respect to the quantization axis parallel 
to the vector r. Integrating over the angle </>' in 
Eq. (|79|l we find that all terms with m' ^ vanish. 
This is because the argument of the Green function 
Ge in Eq. I|79l) depends only on the angle 6'. Using 
the relations Y w (6',0) = i l y/{2l + 1)/4ttP;(cos 0') and 
DoLfa, r , 0) = i~ l y/4n/(2l + l)Yi m (0 r , r ), we then re- 



duce Eq. ifrnji to the form 



I = J2J2 IlmY lm (6 r ,<t> r y, 



(82) 



1=0 m=-l 






/•OO 


Il m = 2ir J d cos 9 


/ r' 2 dr' 




Jo 



/i im (r')P/(cos0') 



-him{r) G E Jr 2 + r' 2 + r\ - 2rr> cosO' . (83) 
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With Eq. I|15(l for the Green function Ge, integration in 
parts over d cos 9' transforms Eq. (|83|) into 



-E r'dr' 



4tt 2 



[hi m (r') - h lm (r)} 



(84) 



The quantities Ai rn represent the sum of integrals in 
which for r_|_ = and r' — > r the integrand either re- 
mains finite or contains a logarithmic integrable singu- 
larity. Setting r+ = one obtains finite values of these 
quantities: 

, _ r r' [h lm (r) + (-l)^h lm (r')}V^E 

•r\-lm — / a o/ i /\ 

J r 47T^{r + r) 

h lm (r')dr' f 1(1 + 1) 



47r 2 r 2 



x Ki (V^E{r + r')) + 
x K (y/=E\r - r'|) + {-l) l K (V^E(r + r') 

- J dxp[' (x)k (y^v 



r 2 + r' 2 — 2rr'x 



(85) 



Equations 1(84(1 and ((85(1 are written for the case E < 0. 
For E > one should put in these equations E instead 
of — E and replace the decaying Bessel functions K\ and 
Kq by the Hankel functions 2iH\ and 2iHo, respectively. 

For calculating the integral in Eq. 1)84(1 we divide the re- 
gion of integration into two parts: inside a small interval 
C, where \r' — r\ <C l/y/\E], and outside this interval. In 
addition, we require that the function hi m (r') does not 
significantly change inside the interval C. The integral 
outside the interval C remains finite and independent of 
r + in the limit r + — > 0. Inside the interval C the argu- 
ment of the Bessel (or Hankel) function is small and the 
integral in Eq. Q84J1 reduces to 



1 

i^ 2 



[hi m (r') - h lm (r)] r' 
r'ec (r' ~ r) 2 + r\ r 



dr'. 



We then expand the quantity \hi m (r') — hi m (r)] in pow- 
ers of (r' — r). Quadratic and higher order terms of the 
expansion give contributions to the integrand, which re- 
main finite for r + — > 0. The contribution of the linear 
term is given by 



1 dhi m {r) 
47r 2 dr 



r'ec 



(r 1 — r) 2 + ' 



-dr'. 



(86) 



For r + —f the integral in Eq. (|86ll remains finite and 
can be written as 



dr' 



r'eC 



(r' — r) 



where the symbol P denotes the principal value of the 
integral. We thus see that in the limit r + — > the whole 



integral / given by Eq. (|79|l is finite. Setting r + = in 
the initial expression (|79|l and keeping the symbol P for 
the integration over dr' , one then writes the integral / in 
the form of Eq. P0|): 

I = P J G E (\r - r'\)[h(r') - h(v)]d\' . 

We now derive Eq. 1)62(1 for the 4-body problem of the 
dimer-dimer scattering and prove that the integral in the 
second line of this equation is convergent. So, changing 
the notation from Y2 to r, we start with the integral 

J = J d 3 r'd 3 R'[h(r',R,')~~ h(r,R)]G(\S - S x \), (87) 

where the argument of the Green function in Eq. (|87|l is 



|5-5 1 | = ^r 2 + (r-r') 2 + (R-R') 2 , 

and consider the limit r\ — > 0. As in the 3-body case, we 
expand the function h(r, R) in spherical harmonics. The 
expansion now reads: 



h(r,K)= E E E hl t'(r,R) 

l—0,l'—0 m— — l m' — — V 



where 8 r , 4> r and #r,, 0r are the polar and azimuthal an- 
gles of the vectors r and R. Then, using Eq. H81|) for 
each of the spherical harmonics in the expansion of the 
function h(r' , R'), we transform Eq. I(87|l to the form 

oo I I 



J = E E E jL n 'Yi m (9rAr)Y lln A0 R ,M:, 



l—O l'—O m— — l m' = — l' 



=4ir 2 f dcos8' r [ dcos9' R I r 
J-i J-i Jo 



(89) 



•OO 

'dr' I R' 2 dR' 
o 



x [h l C (r',i?')fl(cos W(cos^) - h\£ (r,R))G(X), 

(90) 

with the angle between the vectors r' and r denoted as 
9' r , the angle between R' and R as and the argument 
of the Green function written as 

X = yV 2 + r 2 + r 12 - 2rr ' cos 9' r + R 2 + R' 2 - 2 RR> cos 6' R . 

We then separate out the region of integration V, 
where the angles 9' r and 9' R are small, r' is close to r, 
and R' is close to R. The region V is determined by the 
relations: 

\r' - r\ < r < r',r,a; (91) 
\R'-R\ <r <^R',R,a; (92) 
9' T < Vr < min{l, a/VrV}; (93) 
< VR < min{l, a/VRR 1 }, (94) 
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and the quantities ry r , rjn are selected such that they sat- 
isfy the inequalities 



The contribution of the linear terms of the expansion, 
after integrating over R' , can be written as 



rjrV rr' 3> \r' — r\, 
7 ?R V^R 7 > \R'-R\ 



(95) 
(96) 



Outside the region V one can directly put r\ = in the 
argument of the Green function and see that the quanti- 
ties J\™ remain finite. These contributions are denoted 
below as 3\™ . In the region V the argument X of the 
Green function can become equal to n and, in principle, 
this region can be thought of as the one leading to a sin- 
gular behavior of the integral J in the limit r x — > 0. We 
will show that this not the case. 

The Green function G is given by Eq. (|58|l . and in the 
region V the argument of the decaying Bessel function, 
V2X/a, is small. Then, for small angles 9 r and #r, from 
Eq. J5SJ) we obtain 



G- 



15 
32tt 4 



{rl+{r'-r) 2 +(R'-R) 2 ^' r 2 rr'+e'£RR' 



.-7/2 



(97) 



Inside the region V the inequalities 1|93[) and (|94|l allow 
us to put P/(cos6»;) = Pi> (cos 9' R ) = 1 in Eq. (JHOJ- On 
the other hand, owing to the inequalities l|95[) and (|96[) . 
we can put infinity for the upper limits of integration 
over 9' r and 8'^. Then, integrating over the angles inside 
the region V and taking into account the contribution 
j\™ from the configuration space outside this region, 
we reduce Eq. (|§U|) to 



47t 2 rdr 



(r 1 - r) 



r^r' dr' 



un lm 



4n 2 RdR 

To 



: >ev rf + 0' - r) 2 (r 2 + (r' - r) 2 ) 1 / 2 
In 



(r, R) f r'dr' <■ 



r'GV 



r 



TO + V r + ( r ' - r ) 2 



r 



(r 2 + (r 1 - r) 2 )V2 



(99) 



where we put r± = in the second integral as it then 
contains only an integrable logarithmic singularity for 
r' — > r. The first integral on the right hand side of 
Eq. I|99|) also remains finite in the limit of r% — > 0. It 
can be written in the form: 



rodr' 



r'ev (r 2 + (r'-r) 2 y/ 2 
rdr' 



+P 



ro 



r , e y (r'-r) (rg + (r'-r) 2 )i/2 



P 



1 

8^ 



[h i ;™(r>,R')-ht 



r'R' 
,R'ev rR 
'(r,R)} 



dr'dR' 



{r 2 + (r'-r) 2 + {R'-R) 2 f/ 2 ' 



(98) 



For calculating the integral in Eq. (|98|) we expand the 
function h(r' , R') in powers of (r' — r) and (R' — R). 
Quadratic and higher order terms of the expansion lead 
to the integrand which remains finite for r\ — + 0, at least 
after a straightforward integration over one of the vari- 
ables, r' or R'. For example, setting n = quadratic 
terms yield 



with the symbol P denoting the principal value of the 
integral. This means that setting r\ = 0, each of the 
quantities J\™ is equal to the principal value of the in- 
tegral in Eq. (|9()(l . taken with respect to the integration 
over dr' . Note that we could first make the integration 
over dr' and reduce the result to the principal value for 
the integration over dR' . One thus sees that for n — > 
the initial integral J 1)87(1 can be represented in the form 
of Eq. ffEfl: 



d 2 h l f{r,R) 



rgdr' 



8n 2 dr 2 



d 2 til™\r, R) 



8Tr 2 dR 2 
(r,R) 



, £ V [r 2 + {r'-r) 2 y/ 2 
r dR' 

fvev (r 2 + (R'-R) 2 y/ 2 



r) 2 h l ' m ' 



4tt 2 3rd R 



dr'ir' 



r'GV 



rR 



In 



ro + V r o + i r ' - r ) 2 



ro 



(r 2 + (r 1 - r) 



1)1/2 



J = P d 3 r'd 3 R'[h(r' , R') - h{r,R)]G{\S - 5i|), 



where the symbol P denotes the principal value for 
the integration over dr' (or dR'), and \S — Si\ — 

V(r-r') 2 + (R-R') 2 - 
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